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Abstract 

We propose a new "bigravity" model with two positive tension AdS^ branes in AdS^ 
bulk and no negative tension branes. The bounce of the "warp" factor mimics the effect 
of a negative brane and thus gives rise to an anomalously light graviton KK mode. This 
configuration satisfies the weak energy condition and has no ghost state. In addition, 
the extra polarization states of the massive graviton practically decouple and thus it does 
not contradict to Einsteinian gravity. However, the model has certain phenomenological 
difficulties associated with the presence of a negative cosmological constant on the branes. 
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1 Introduction 



Brane universe models in more than four dimensions have been extensively studied 
over the last three years because they provide novel ideas for resolutions of long standing 
problems to particle physics such as the Planck hierarchy ^ one. Moreover, mechanisms 
to localize gravity on a brane ^ have led to the realization that if extra dimensions 
exist, they need not be compact. Also "bulk" (transverse to the 3-brane space) physics 
turns out to be very interesting giving alternative possible explanations to other puzzles 
of particle physics, like the smallness of the neutrino masses, the neutrino oscillations, or 
the pattern of the SM fermion mass hierarchy (see for example and references therein). 
The phenomenological implication of these constructions are radical and the fact that these 
ideas can, in principle, be put to test in current and future experiments, makes them very 
attractive. A comprehensive account of these ideas can be found in 0]. 

In this paper, we present a model which belongs to a class of brane universe models 
P, ^ |lT| that suggest that a part or all of gravitational interactions come from massive 
gravitons. In the first case ( "bigravity" ) , gravitational interactions are the net effect of 
a massless graviton and a finite number of KK states that have sufficiently small mass 
(or /and coupling) so that there is no conflict with phenomenology. In the second case 
( "multigravity" ) , in the absence of a massless graviton, the normal Newtonian gravity at 
intermediate scales is reproduced by special properties of the lowest part of the KK tower 
(with mass close to zero). These models predict that at sufficiently large scales, which 
correspond to the Compton wavelength of the KK states involved, modifications to either 
the Newtonian coupling constant or the inverse square law will appear due to the Yukawa 
suppression of the KK states contribution. These modifications affect the CMB power 



spectrum and thus are, in principle, testable |12 



The prototype model of this class was the " H h" "bigravity" model ^ where we 

considered a modification of the compact Randall- Sundrum model (RSI) with two pos- 
itive tension fiat branes ("-I-" branes) separated by one intermediate negative tension fiat 
brane ("— " brane) in AdS^ bulk. The task of finding the KK spectrum reduces to a simple 
quantum mechanical problem. It is simple to see that the model (as every compact model) 
has a massless graviton that corresponds to the ground state of the system whose wave- 
function follows the "warp" factor. It is also easy to see as in figure 1 that there should be a 
state with wavefunction antisymmetric with respect to the minimum of the "warp" factor, 
whose mass splitting from the massless graviton will be very small compared to the ones 
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Figure 1: The graviton (solid line), first (dashed line) and second (dotted line) KK states 

wavefunctions in the symmetric " H h" model. The wavefunctions are not smooth on the 

"— " branes. The same pattern have also the " + +" model wavefunctions with the position 
of the " brane corresponding to the minimum of the warp factor. The wavefunctions are 
then smooth. 



of the higher levels. Because the "warp" factor is exponential the difference of the mass 
behaviour of the first and the rest of the KK states is also exponential. This allowed for 
the construction of a "bigravity" model in which the remainder of the KK tower does not 
affect gravity beyond the millimeter bound. In this case gravity at large scales beyond the 
Compton wavelength of the first state will be modified in the sense that Newton's constant 
will be reduced. In particular, it is reduced by one half in the symmetric configuration and 
it can be effectively switched off in the highly asymmetric case. 

Independently, Gregory, Rubakov and Sibiryakov (GRS) suggested a construction in 
which gravity is also modified at ultra-large scales. The GRS model modifies the decom- 
pactified Randall-Sundrum model (RS2) by adding a "— " brane of half the tension of the 
brane and requiring flat space to the right of the new " brane. This model does not 
have a normalizable 4D graviton but generates 4D gravity at intermediate distances due to 



a resonance-like behaviour |T3|, of the wavefunctions of the KK states continuum. Grav- 
ity in this picture appears to be "quasi-localized" and this time it is modified at ultra-large 
scales in the sense that the inverse square law becomes inverse cube, i.e. five dimensional. 



Although the " H h" and GRS models look quite different it was shown in |[T^ and [|TT 

that these two models are the limiting cases of a more general " H h" multi-brane model 

that interpolates between the "bigravity" " H h" model and the "quasi-localized" gravity 
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GRS model. The key element of both models that made "multigravity" possible was the 
bounce of the "warp" factor generated by the "— " branes. 

However, one should be careful when dealing with models with ultralight massive KK 
states because it is known that the extra polarizations of the massive gravitons do not 



decouple in the limit of vanishing mass, the famous van Dam - Veltman - Zakharov [15 



discontinuity. This could make these models disagree [1^1 with standard tests of General 
Relativity, as for example the bending of the light by the sun. Furthermore, the moduli 
(radions) associated with the perturbations of the "— " branes are necessarily physical ghost 
fields |]16|, therefore unacceptable. The latter problem is connected to the violation of the 
weaker energy condition on " branes sandwiched between branes. In the GRS 



model this radion cancels |T8| the extra polarizations of the massive gravitons and gives the 
graviton propagator the correct tensorial structure at intermediate distances. However, the 
model has still an explicit ghost in the spectrum which reveals itself as scalar antigravity 
at cosmological scales. A mechanism of cancelling both the extra massive graviton polar- 



izations and the radion field contribution was suggested in |10, 11 1 and involves some bulk 



dynamics which are necessary to stabilize the system, based on a scenario described in 
This mechanism is however non-local in the extra dimension and because of this may not 
be very attractive. 

In the present paper we will demonstrate that there is actually a way out of both these 
problems. We will use a two brane model with only "+" branes which was known in 
|2T|, |2^, ^ to exhibit a bounce of the "warp" factor, and therefore is bound 

to have "bigravity" by general arguments presented in ||^, This can be achieved if we 
consider two AdS^ branes in AdS^ bulk. Motivated by quasi-localization of gravity, Randall 
and Karch studied the non-compact case of a single "+" brane, which is a limiting case 
of the asymmetric " + +" model. The weaker energy condition is satisfied so there is no ghost 
modulus in this setup. Furthermore, as was shown in in AdS space it is possible to 

circumvent the van Dam - Veltman - Zakharov no go theorem about the non-decoupling of 
the massive graviton extra polarization states. The price we pay is that there is a remnant 
negative cosmological constant on the brane. This sets an horizon scale and unfortunately 
the Compton wavelength of the light state of the system lies exponentially far from this 
horizon. This does not change even if we consider a highly asymmetric version of this 
model. As a result, this "bigravity" model makes no predictions for observable deviations 



''One day after this paper had appeared in the hep-archives, ref.[p9[ appeared, reaching the same con- 
clusions by a different method. 
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from Newtonian gravity at ultra-large distances. In addition, although theoretically there 
exist modifications of General Relativity at all scales due to the additional polarization 
states of the massive graviton, they are so highly suppressed that they are not observable. 
Additionally, let us note that recently models with only positive tension branes and similar 
"warp" factors (created in a different and dynamical way) were discussed in [^, E^, |3^ . 



2 The two positive brane model 

The model consists of two 3-branes with tensions Vi and V2 respectively, in an AdS^ space 
with five dimensional cosmological constant A < 0. The 5-th dimension has the geometry 
of an orbifold and the branes are located at its fixed points, i.e. Lq = and Li = L. Due 
to orbifolding, we can restrict ourselves to the region < z < L, imposing the suitable 
boundary conditions to our solutions. Firstly, we find a suitable vacuum solution. The 
action of this setup is: 

S= [ d^x r rfzv^{-A + 2M'^R} + d^xV^^-Gi^ (1) 

J J—L ■ J z=Li 

where G^^l is the induced metric on the branes. The notation is the same as in Ref. 
The Einstein equations that arise from this action are: 

\ 

) 

/ 



1 1 ( y— 

RmN - -GmnR = I ^^MN - X! ^JZlQ ^'t^l^M^N^i^ ~ ^i) 



(2) 



In order to find a specific form of the equations of motion we need to write a metric 
ansatz which will take in account the spacetime symmetries of the 3-brane. Since we would 
like not to restrict our model to fiat solutions on the branes, we should make a choice which 
will let us interpolate between the maximally symmetric space-times in four dimensions. 



i.e. de-Sitter, Minkowski and Anti-de-Sitter. The metric ansatz |jT^ that accomplishes this 
is the following: 

ds^ = a\z){-de + e'^^'dx^) + h\z)dz'' (3) 

where H is the "Hubble" parameter and is determined in terms of the brane tension Vi and 
the bulk cosmological constant A from Einstein's equations. The z-dependent function a{z) 
is the "warp" factor that is essential for the localization of gravity and also for producing 
the hierarchy between the two branes. In the case of fiat brane solution, i.e. the effective 
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cosmological constant on the brane is zero, we have H — 0. On the other hand if we 
demand a de-Sitter solution on the brane, i.e. the effective cosmological constant on the 
branes is positive, we have H"^ > 0. In the case of Anti-de-Sitter solution, i.e. the effective 
cosmological constant on the branes is negative, we have < and thus H is imaginary. 
In order to get a physical interpretation of the latter case it is necessary to analytically 
continue the solution by a coordinate transformation of the form t — —ix'i, xi — it', 
X2 — x'2 and x-i — x'^. After this transformation the metric ansatz can be written in the 
following form: 

ds" = a'{z){dx'l + e^^^ + ^^^2 ^ ^^/^2^^ ^ \?{z)dz^ (4) 

Furthermore, in order to have a more compact notation for all cases of maximally 
symmetric spaces and simplify our calculations, it is useful to bring the metric ansatz in 
the form: 



ds = TTT^'n^^dx^dx'' + h {z)dz (5) 

where = rj^^x^x^ . It can be shown that the Ricci scalar for this metric is i? = —12H^. 
Thus this metric represents all maximally symmetric spaces: Minkowski for = 0, Anti- 
de-Sitter for > and de-Sitter for < 0. From now on we shall choose the gauge 
b{z) = 1 where our coordinate system is Gaussian Normal. A straightforward calculation 
of the Einstein's equations gives us the following differential equations for a{z): 

i 

By solving the above equations we find that the solution can be written in the form : 

Vik 

a{z) — cosh(A;|2;|) -| — — sinh(/i;|2;|) (8) 



A 



with 



j^{V^^k^-A^) , ^ < 1/1 for (i54 branes 

for flat branes (9) 

g(A2 - V^e) , ^ > 14 for AdS^ branes 
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where we have normahzed a(0) = 1 and assumed Vi > 0. Also we have defined k = \J -^aW- 

Additionally, in order to have this solution, the brane tensions Vi, V2, the bulk cosmo- 
logical constant |A| and the position of the second brane L must be related through the 
equation: 

Let us now restrict ourselves to the case of AdS/^ spacetime on the two branes which will 
turn out to be the most interesting. In this case the condition ^ > Vi must hold. Hence, 
we can define tanh(fczo) = ^ and write the solution in the form: 

a{z) = ^^^MM^o-N)) (11) 
cosh(K2;o) 

from which it is clear that the "warp" factor has a minimum aX z = Zq. From this point 
we can see the role of the 74^5*4 on the branes, i.e. the role of the condition ^ > Vi. This 
condition allows us to have the bounce form of the "warp" factor {i.e. a minimum in the 
"warp" factor) allowing the second brane to have positive tension and give us, as we will 

see shortly, a phenomenology quite similar to the " H V" "bigravity" model This can 

be easily seen from the eq.([TO|) which relates the brane tensions and the distance between 
the branes. From this equation we indeed see that by placing the second brane after the 
minimum of the "warp" factor we can make the tension of the second brane positive and 
thus both branes that appear in the model have positive tension avoiding the problems 
associated with negative tension branes. In fact it is clear that the present model mimics 

the characteristics of the " H h" "bigravity" model since what we effectively do is to 

reproduce the effect of the presence of a negative tension brane, i.e. the bounce form of the 
"warp" factor, with another mechanism allowing a negative four dimensional cosmological 

constant on the brane. Note that in the limit that ^1 — * ^ (fiat limit) the minimum of 
the "warp" factor tends to infinity and if we wish to have a brane at a finite point, it will 
necessarily have negative tension. 

The relationship between the 4D effective fundamental scale and the five dimen- 
sional fundamental scale M can be easily found by dimensional reduction to be: 

^the factor 2M^ multiplies the four dimensional Ricci scalar in the Lagrangian after dimensionally 
reducing 
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= 27 7 [kL + smh{kL) cosh(A;(L - 2^0))] (12) 

/ccosh {kzQ) 

The above formula tells us that for finite L the compactification volume is finite and 
thus the zero mode is normalizable. In the case where we send the second brane to infinity, 
the compactification volume becomes infinite which indicates that the zero mode becomes 
non-normalizable. Note that is not necessarily equal to Mpi since as will see shortly, at 
least for a sector of the parameter space of our model, gravity is the result not only of the 
massless graviton but also of an ultralight KK state. 

The "warp" factor renormalizes the physical scales of the theory as in [^]. Thus, all 
mass parameters mo on a brane placed at the point z are rescaled as 



m = a{z)mQ (13) 

Hence, assuming some kind of stabilization mechanism which fixes the positions of the 
branes, one can choose a distance between the two branes such that this rescaling leads to 
the creation of a desired mass hierarchy. 

However, since we consider non-fiat solutions on the branes, we have to make sure that 
the four dimensional effective cosmological constant does not contradict present experi- 
mental and observational bounds. Recent experimental data favour a positive cosmological 
constant, nevertheless since zero cosmological constant is not completely ruled out it can 
be argued that also a tiny negative cosmological constant can be acceptable within the 
experimental uncertainties. The effective cosmological constant on the two branes is: 
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Am = -l2HHdl = ^-^^eul (14) 

cosh {kzQ) 

From the above formula we can see that we can make the cosmological constant small 
enough |A4rf| < lO^^^^Mpj if we choose large enough kz^, i.e. kzo > 135. This however will 
make observable deviations from Newtonian gravity at ultra-large scales impossible as we 
will see in the next section. 

To determine the phenomenology of the model we need to know the KK spectrum that 
follows from the dimensional reduction. This is determined by considering the (linear) 
fluctuations of the metric around the vacuum solution that we found above. We can write 
the metric perturbation in the form: 
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(ix'^dx" + dz'^ (15) 



where g^^f^ is the vacuum solution. Here we have ignored the radion mode that could 
be used to stabilize the brane positions z — Lq and z — Li, assuming some stabilization 
mechanism. We expand the field h^^{x, z) into graviton and KK plane waves: 

oo 

V(^,^) = E/^£H^)^^"n^) (16) 

n=0 

where we demand (V^V" + 2H'^ — rri^) h^^^J = and additionally V^h^^^ = h^"l°' = 0. The 

function "^^"'\z) will obey a second order differential equation which after a change of 
variables and a redefinition of the wavefunction reduces to an ordinary Schrodinger-type 
equation: 

{—91 + V{w)] ¥-\w) = '^¥-\w) (17) 
where the potential is given by: 



V{w) = - — + 



8 cos^ (^(I'W^I ~ ""^o) 



3k 
T 



cosh.{kzo) 



(18) 



with k defined as k 



cosh(fczo) 

are related with the old ones by: 



The new variables and the redefinition of the wavefunction 



w = sgn(2;) ~ 
k 



arctan [ tanh(— ^LL — —)] + arctan [ tanh(^^^ 



\/a{zj 



(19) 



(20) 



Thus in terms of the new coordinates, the branes arc at = and w^, with the mini- 
mum of the potential at Wq = r arctan (tanh(^)). Also note that with this transformation 



the point 2; = oo is mapped to the finite point 
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_ 2 
°° k 



f + arctan (tanh(^)) 



From now on we restrict ourselves to the sjTumetric configuration of tlie two branes 
witli respect to tfie minimum wo [i.e. the first brane at and the second at 2wo ), since 
the important characteristics of the model appear independently of the details of the con- 
figuration. Thus, the model has been reduced to a "quantum mechanical problem" with 
(5-function potentials wells of the same weight and an extra smoothing term in-between 
(due to the AdS geometry). This gives the potential a double "volcano" form. 

An interesting characteristic of this potential is that it always (for the compact cases 
i.e. wl < Woo) gives rise to a normalizable massless zero mode, something that is expected 
since the volume of the extra dimension is finite. The zero mode wavefunction is given by: 

¥'^\w) = . — (21) 

[cos(A;(wo — |w|))]'^/^ 



dw '^^^'(w) = 1 

-WL 

chosen so that we get the standard form of the Fierz-Pauli Lagrangian. 

The form of the zero mode resembles the one of the zero mode of the " H h" model, 

i.e. it has a bounce form with the turning at Wq (see figure 1). In the case of the " H h" 

the cause for this was the presence of the "— " brane. In the present model it turns out that 
by considering AdS spacetime on the branes we get the same effect. 

In the case that we send the second brane to infinity {i.e. w —>■ Woo) the zero mode 
fails to be normalizable due to singularity of the wavefunction exactly at that point. This 



w 



can be also seen from eq.(|T^) which implies that at this limit M* becomes infinite {i.e. the 
coupling of the zero mode becomes zero). Thus in this limit the model has no zero mode 
and all gravitational interactions must be produced by the ultralight first KK mode. The 



spectrum in this case was discussed by Randall and Karch at . 

Considering the Schrodinger equation for m 7^ we can determine the wavefunctions of 
the KK tower. It turns out that the differential equation can be brought to a hypergeometric 
form, and hence the general solution is given in terms two hypergeometric functions: 



cos 



'/^{k{\w\ - Wo)) \C, F{d^, K, i; sm\~k{\w\ - Wo))) 



+ C2 I sin(A;(|ti?| - Wq))\ F{an + |, &n + |, |; sin^(A;(|w| 



^0))) 



(22) 
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where 




The boundary conditions (i.e. the jump of the wave function at the points w — 0, wl) 
result in a 2 X 2 homogeneous hnear system which, in order to have a non-trivial solution, 
leads to the vanishing determinant. In the symmetric configuration which we consider, 
this procedure can be simplified by considering even and odd functions with respect to the 
minimum of the potential Wq. 

In more detail, the odd eigenfunctions obeying the b.c. '^^'^\wq) = will have Ci = 
and thus the form: 

^W = C2COS^/2^^(|«;|-«;o))|sin(fc(|^/;|-^/;o))|F(5, + ^,6, + i,^;sin^ (24) 

On the other hand, the even eigenfunctions obeying the b.c. '^^"'^ 'i'^o) = will have 
C2 — and thus the form: 



^CiCOs'/''{^\w\-Wo))F{d„X,\;sm\^\w\-Wo))) (25) 
The remaining boundary condition is given by: 

'(0) + — tanh(A;zo)*^''^ (0) = (26) 

and determines the mass spectrum of the KK states. From this quantization condition 

we get that the KK spectrum has a special first mode similar to the one of the " H V" 

"bigravity" model. For kz^ ^ 5 the mass of the first mode is given by the approximate 
relation: 

mi = 4^3 A; e'^*"" (27) 

In contrast, the masses of the next levels, if we put together the results for even and 
odd wavef unctions, are given by the formula: 

m„+i = 2yn(nT3) k e'^^^ (28) 
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with n = 1,2, .... 

We note that the first KK state has a different scaling law with respect to the position of 
the minimum of the "warp" factor compared to the rest of the KK tower, since it scales as 
g~2fc^o -while the rest of the tower scales as e~^^°. Thus the first KK state is generally much 
lighter than the rest of the tower. It is clear that this mass spectrum resembles the one of 

the "H h" "bigravity" model. The deeper reason for this is again the common form of the 

"warp" factor. In both cases the "warp" factor has a minimum due to its "bounce" form. 
The graviton wave function follows the form of the "warp" factor, i.e. it is symmetric with 
respect to Wq, while the wavefunction of the first KK state is antisymmetric in respect to 
Wo (see figure 1). The absolute values of the two wavefunctions are almost identical in all 
regions except near wq (the symmetric is nonzero while the antisymmetric is zero at Wq). 
The graviton wavefunction is suppressed by the factor ^^^.^l^^j^^^-^ at Wq which brings it's value 

close to zero for reasonable values of kzQ. Thus, the mass difference which is determined 
by the wavefunction near wq is expected to be generally very small, a fact which formally 
appears as the extra suppression factor e~^^° in the formula of mi in comparison with the 
rest of the KK tower. 

In the case that we consider an asymmetric brane configuration, for example wl> 2wo 
the spectrum is effectively independent of the position of the second brane wl (considering 
kzQ > 5). Thus, even in the case that we place the second brane at Woo, i.e. the point which 
corresponds to infinity in the 2;-coordinates, the spectrum is given essentially by the same 
formulas. In the case that the second brane is placed at wq < wl < ^wq, some dependence 
on the position of the second brane {i.e. dependence on the scale hierarchy between the 
branes) is present. Nevertheless, the main characteristics of the spectrum remain the same, 
i.e. the first KK state is special and always much lighter than the others. In conclusion, 
the key parameter which determines the spectrum is the position of the minimum of the 
"warp" factor. 

Returning to our wavefunction solutions, we should note that for each eigenfunction 
the normalization constants Ci and C2 can be determined by the normalization condition 

/ dw "^^^'{w) = 1 which is such that we get the standard form of the Fierz-Pauli 

Lagrangian for the KK states. Knowing the normalization of the wavefunctions, it is 
straightforward to calculate the strength of the interaction of the KK states with the SM 
fields confined on the braneQ. This can be calculated by expanding the minimal gravitational 
^In the symmetric configuration it does not make any difference wliicti brane is our universe. 
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couplmg of the SM Lagxa„g,au / d'x (6. SM f,Ms) wth respect to the metr.e. In 

this way we get: 

n>0 

= -J^,h^»T,. (X) - E '''"'^!;r^ enx)T,. (X) (29) 

with T^i, the energy momentum tensor of the SM Lagrangian. Thus, the couphng of the 
zero and the first KK mode to matter are respectively: 



1 A 1 



Co M3/2 

Cl M3/2 ~ M, 



(30) 
(31) 



where A is the zero mode normahzation constant which turns out to be "^^fj— • We should 
also note that the couplings of the rest of the KK states are much smaller and scale as 

Exploiting the different mass scaling of the first KK relative to the rest we can ask 

whether it is possible to realize a "bigravity" scenario similar to that in " H h" "bigravity" 

model. In that model by appropriately choosing the position of the minimum of the "warp" 
factor, it was possible to make the first KK state have mass such that the corresponding 
wavelength is of the order of the cosmological scale that gravity has been tested and at 
the same time have the rest of the KK tower wavelengths below 1mm (so that there is no 
conflict with Cavendish bounds). In this scenario the gravitational interactions are due to 
the net effect of the massless graviton and the first ultralight KK state. From eq.(^), (pi|) 
it can be understood that in the symmetric configuration the massless graviton and the 
special KK state contribute by the same amount to the gravitational interactions. In other 
words: 

' ^ +4 - M„^^| (32) 



In the present model, the fact that the effective four dimensional cosmological constant 
should be set very close to zero, requires that the "warp" factor is constrained by kz^ > 135 
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and thus, in this case, the spectrum of the KK states will be very dense (tending to con- 
tinuum) bringing more states close to zero mass. The KK states that have masses which 
correspond to wavelengths larger than 1mm have sufficiently small coupling so that there is 
no conflict with phenomenology (the situation is exactly similar to the RS2 case where the 
coupling of the KK states is proportional to their mass and thus it is decreasing for lighter 
KK states). The fact that the spectrum tends to a continuum shadows the special role of 
the ffist KK state. Moreover, it is interesting to note that at the limit where the minimum 
of the "warp" factor is sent to infinity (woo) the special behaviour of the ffist KK persists 



and does not catch the other levels (by changing its scaling law) as was the case in 
This means that the limit w Woo will indeed be identical to two RS2, but on the other 
hand it is interesting to note that what we call graviton in the RS2 limit is actually the 
combination of a massless graviton and the "massless" limit of the special first KK state. 
This "massless" limit exists as we will see in the next section and ensures that locality is 
respected by the model, since physics on the brane does not get affected from the boundary 
condition at infinity. 

3 Discussion and conclusions 

The fact that we have a ultralight graviton in our spectrum is at first sight worrying because 
it is well known that in the flat space the tensor structure of the propagators of the massless 
and of the massive graviton are different [113 1 and that there is no smooth limit between 



them when m — 0. The bending of the light by the sun agrees with the prediction of 
the Einstein theory to 1% accuracy. This is sufficient to rule out any scenario which a 
significant component of gravity is due to a massive graviton, however light its mass could 



be. However, as was shown in [|28[, the situation in AdS space is quite different. There it 
was shown that if we could arrange ^ < 0.1 there is no discrepancy with standard tests 
of Einsteinian gravity as for example the bending of the light by the sun. 

In the particular model we have at hand, it is ^ ~ g-^^o gQ gg^j^ easily accommodate 
the above bound. Then, the Euclidean propagator (in configuration space) of the massive 
KK states for relatively large zq will be given by: 

G;r.;^...(x, y) = ^^('^MM"^--' + W - (l - ^e-''^"") 5^J^.,.) (33) 
where /x is the geodesic distance between two points. In the above, we have omitted terms 
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that do not contribute when integrated with a conserved T^y. For kz^ > 2.3 there in no 
problem with the bending of hght. However, if our aim is to see modifications of gravity at 
ultra-large distances, this is impossible because the Compton wavelength of our ultralight 
graviton will be e'^^" times bigger than the horizon of the AdS^^ space on our brane due 
to equations (|27D , (pSf ). The "Hubble" parameter follows m2 rather than mi. 

What happens if one takes an asymmetric version of this model where L > 2zq is that 
the spectrum does not get significantly modified so we are effectively in the same situation. 

In the case where zq < L < 2zq the spectrum will behave like the one of the " H h" model. 

Then for ^ 1 we will have mi ~ tue'^^^M ~ e'^'^'Mpi, H ~ e~^^°M ~ e"''^Mpi where 
M «i M-pi/uj is the fundamental scale, uj ~ e-(2^o-i) jg ^j-^g "warp" factor and x = L — zq. 
Again, the ultralight graviton is hiding well beyond the AdS horizon. However, the coupling 
of the remaining of the KK tower to matter will be different than the symmetric case and 
one may have different corrections to Newton's law on the left and right branes. 

In summary, in this paper we presented a new "bigravity" model with two AdS^ branes 

in AdS^ bulk which has a lot of similarities with the " H h" "bigravity" model. The fact 

that we have no "— " branes removes the ghost state problem and furthermore, due to some 
amazing property of the AdS space we are able to circumvent the van Dam - Veltman - 
Zakharov discontinuity of the graviton propagator. This makes the model compatible with 
the predictions of General Relativity in the small graviton mass limit since the extra degrees 
of freedom of the massive graviton practically decouple. However, the presence of the AdS 
horizon prevents the modifications of gravity at large distances becoming observable. In the 
future it would be interesting to explore similar models where it would be possible to obtain 
observable modifications of gravity at large (cosmological) scales. In particular, it would 
be interesting to see if the above characteristics of this model persist when we add matter 
density on the branes. Then it is interesting to examine if /how the ultralight graviton is 
going to reveal itself in the cosmological solutions discussed by ||25|| . 
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Addendum: One day after this work had appeared in the hep-archives, ref. 
appeared, being the pubhshed version of In that paper the non-compact case of a 



single brane was studied. 



References 



[1] K. Akama in Gauge Theory and Gravitation, Proceedings of the International Sympo- 
sium, Nara, Japan, 1982, ed. K.Kikkawa, N.Nakanishi and H. Nariai (Springer- Verlag, 
1983), 267; e-version: K.Akama, |Eep-th/0001lT^ ; 



V.A.Rubakov and M.E. Shaposhnikov, Phys. Lett. B125 (1983) 136; 
M. Visser, Phys. Lett. B159 (1985) 22; 
E. J. Squires, Phys. Lett. B167 (1985) 286. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. Dvah, Phys. Lett. B429 (1998) 263; 
N. Arkani-Hamed, S. Dimopoulos and G. Dvah, Phys. Rev. D59 (1999) 086004; 
I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 (1998) 
257. 

[3] M. Gogberashvih, Mod. Phys. Lett. A14 (1999) 2025. 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370. 

[5] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690. 



[6] N. Arkani-Hamed, S. Dimopoulos, G. Dvali and J. March-Russell, Piep-ph /98 1 1448 
Y. Grossman and M. Neubert, Phys. Lett. B474 (2000) 361; 
G. Dvali and M. Shifman, Phys. Lett. B475 (2000) 295. 

[7] Online proceedings of the conference "Physics beyond four dimensions", July 2000, 
ICTP, Trieste, www.ictp.trieste.it/~highener/beyond4d.html . 

[8] I.I. Kogan, S. Mouslopoulos, A. Papazoglou, G.G. Ross and J. Santiago, Nucl. Phys. 
B584 (2000) 313; 

S. Mouslopoulos and A. Papazoglou, JHEP 0011 (2000) 018. 



15 



[9] R. Gregory, V.A. Rubakov and S.M. Sibiryakov, Phys. Rev. Lett. 84 (2000) 5928. 
[10] I.I. Kogan and G.G. Ross, Phys. Lett. B485 (2000) 255. 



[11] I.I. Kogan, S. Mouslopoulos, A. Papazoglou and G.G. Ross, [hep-th/0006030| , to appear 
in Nucl. Phys. B. 



[12] P. Binetruy and J. Silk, |astro-ph/ 0007452 



[13] C. Csaki, J. Erlich and T.J. Hollowood, Phys. Rev. Lett. 84 (2000) 5932. 

[14] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B484 (2000) 112. 

[15] H. van Dam and M. Veltman, Nucl. Phys. B22 (1970) 397; 
V.I. Zakharov, JETP Lett. 12 (1970) 312. 

[16] G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B484 (2000) 129; 
L. Pilo, R. Rattazzi and A. Zaffaroni, JHEP 0007 (2000) 056. 



[17] D.Z. Freedman, S.S. Gubser, K. Pilch and N.P. Warner, |hep-th/99040T7 



E. Witten, |hep-ph/ 0002297 . 



C. Csaki, J. Erlich and T.J. Hollowood, Phys. Lett. B481 (2000) 107; 
R. Gregory, V.A. Rubakov and S.M. Sibiryakov, Phys. Lett. B489 (2000) 203. 

[19] N. Kaloper, Phys. Rev. D60 (1999) 123506. 

[20] T. Nihei, Phys. Lett. B465 (1999) 81. 

[21] H.B. Kim and H.D. Kim, Phys. Rev. D61 (2000) 064003. 

[22] O. DeWolfe, D.Z. Freedman, S.S. Gubser and A. Karch, Phys. Rev. D62 (2000) 046008. 

[23] P. Kanti, I.I. Kogan, K.A. Olive and M. Pospelov, Phys. Rev. D61 (2000) 106004. 

[24] P. Kanti, K.A. Olive and M. Pospelov, Phys. Lett. 481 (2000) 386. 

[25] P. Kanti, K.A. Olive and M. Pospelov, Phys. Rev. D62 (2000) 126004. 

[26] Talk by I.I. Kogan, www.ictp.trieste.it/~highener/beyond4d.html . 



16 



[27] Talks by L. Randall, www.ictp.trieste.it/~liigliener/beyond4d.html ; 
and A. Karch, liep.pliysics.uoc.gr/string2000/programl.html . 



[28] I.I. Kogan, S. Mouslopoulos and A. Papazoglou, |hep-th/0011138 
[29] M. Porrati, |F^-th/0011152 . 

[30] John EUis, N.E. Mavromatos, D.V. Nanopoulos, Phys. Rev. D62 (2000) 084019; 
G.K. Leontaris and N.E. Mavromatos, Phys. Rev. D61 (2000) 124004; 
G.K. Leontaris and N.E. Mavromatos, |hep-th/OOlTTO^ . 

[31] S. Nojiri and S.D. Odintsov, Phys. Lett. B484 (2000) 119; 
S. Nojiri, S.D. Odintsov and K.E. Osetrin, |hep-th/ 00090591 . 

[32] A.S. Gorsky and K.G. Selivanov, Phys. Lett. B485 (2000) 271; 



A.S. Gorsky and K.G. Selivanov, |hep-th/0006044 



A. Karch and L. Randall, |hep-th/001lT56 



17 



